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We probe the time dependence of the dark energy equation of state (EOS) in light of three-year 
WMAP (WMAP3) and the combination with other tentative cosmological observations from galaxy 
clustering (SDSS) and Type la Supernova (SNIa). We mainly focus on cases where the EOS is 
oscillating or with local bumps. By performing a global analysis with the Markov Chain Monte 
Carlo (MCMC) method, we find the current observations, in particular the WMAP3 + SDSS data 
combination, allow large oscillations of the EOS which can leave oscillating features on the (residual) 
Hubble diagram, and such oscillations are potentially detectable by future observations like SNAP, 
or even by the CURRENTLY ONGOING SNIa observations. Local bumps of dark energy EOS can 
also leave imprints on CMB, LSS and SNIa. In cases where the bumps take place at low redshifts and 
the effective EOS is close to —1, CMB and LSS observations cannot give stringent constraints on such 
possibilities. However, geometrical observations like (future) SNIa can possibly detect such features. 
On the other hand when the local bumps take place at higher redshifts beyond the detectability 
of SNIa, future precise observations like Gamma-ray bursts, CMB and LSS may possibly detect 
such features. In particular, we find that bump-like dark energy EOS on high redshifts might be 
responsible for the localized features of WMAP on ranges I ~ 20 — 40, which is interesting and 
deserves addressing further. 



PACS number(s): 98.80.Es, 98.80.Cq 



I. INTRODUCTION 

The three year Wilkinson Microwave Anisotropy Probe observations (WMAP3)Q, 0, 0, 0, HI have made so far 
the most precise probe on the Cosmic Microwave Background (CMB) Radiations. In the fittings to a constant 
equation of state (EOS) of dark energy (DE) w, combinations of WMAP with other cosmological observations are 
in remarkable agreement with a cosmological constant (CC) except for the WMAP + SDSS combination, where 
w > —1 is favored a bit more than lcrQ. The measurements of the SDSS power spectrum[(| Q in some sense 
make the most precise probe of the current linear galaxy matter power spectrum and will hopefully get significantly 
improved within the coming few years. If the preference of w > — 1 holds on with the accumulation of cosmological 
observations this will also help significantly on our understandings towards dark energy. A cosmological constant, 
which is theoretically problematic at present 8, 9], will NOT be the source driving the current accelerated expansion 
and a preferred candidate would be something like quintessence jlOj. On the other hand, the observations from the 
Type la Supernova (SNIa) in some sense make the only direct detection of dark energy^l 0, 0, 0, 0, and 
currently a combination of WMAP + SNIa or CMB + SNIa + LSS are well consistent with the cosmological constant 
and the preference of a quintessence-like equation state has disappeared Q. It is noteworthy that in the combinations 
with the Lyman a forest Ref.0] shows that a constant EOS w < —1 is preferred slightly. Moreover when one 
considers the observational imprints by dynamical equation of state, an EOS which gets across —1 is mildly favored 
by the current observations [13 • Intriguingly, we are also aware that the predictions for the luminosity distance- 
redshift relationship from the ACDM model by WMAP only are in notable discrepancies with the "gold" samples 
reported by Riess et aZ|l4| . Although the discrepancy might be due to some systematical uncertainties in the Riess 
"gold" sample^H, this needs to be confronted with the accumulation of the 5-year SNLS observations and the 
ongoing SNIa projects like the Supernova Cosmology Project (SCP) and from the Supernova Search Team (SST). 
Alternatively, this might be due to the implications of dynamical dark energy with oscillating equation of state. 

Although the temperature-temperature correlation (TT) power of WMAP3 is now cosmic variance limite d up to 
I ~ 400 and the third peak is now detected, the tentative features as discovered by the first year WMAP |20l l2ll |22| 
are still present: the low TT quadrupole and localized oscillating features on TT for I ~ 30 - 500- Although the 
signatures of glitches on the first peak as discovered by the first year WMAP have now become weak, they do exist 
and go beyond the limited cosmic variance [a. While for the low WMAP TT quadrupole many authors are inclined 
to attribute it to cutoff primordial spectrum|23l |24| , and even BEFORE the release of the first year WMAP Ref. [2{| 
claimed oscillating primordial spectrum could lead to oscillations around the first peak of CMB TT power, similar 
effects might be due to features on dark energy rather than inflation. For example, Ref. 26] has attributed the low 
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quadrupole to some subtle physics of dark energy during inflation. 

In the literature there have been many investigations on inflationary models with broken scale invariance|29l |30L [3 



WMAP data could not exclude a large running of the spectral index (for relevant study see also [43), which has been 
somewhat dramatically confirmed by the first year WMAP and WMAP3 in combination with other observations 0,0 
except for the case with the Lyman alpha forest 0, Inflation and dark energy, both of which describe the 

accelerated ex pan sion of the universe, might have some relations and Ref. 48] proposed a new picture of quintessential 
inflation. Ref. 49] has made an attempt trying to find such relations from the observational aspect. While oscillating 
primordial spectrum may be responsible for the glitches on CMB, oscillating EOS of dark energy may be helpful to 
solve the coincidence problem of dark energy[S(],[in],[52,[5i|. In Ref.p5l| in the framework of QuintomjS^, an attempt 
was carried out to unify dark energy and inflation, meanwhile solving the coincidence problem of DE. 

On theoretical aspect dark energy is among the biggest problem of modern cosmology [si Wi IHH l56l l57l l5St . Dy- 
namical dark energy models rather than the simple cosmological constant have attracted more interests in theoretical 
studies 6(3] • In cases where the mysterious component of DE is driven by scalar fields0,|H3,IH3 the EOS is typically 
not like that by a cosmological constant, this opens a possibility for us to tell CC from scalar dark energy models with 
the cosmological observations. Moreover the current observations have already opened a robust window to probe the 
behavior of dark energy independently, and in cases when w ^ —1 is preferred, dynamical DE models which satisfy 
the observations are put forward 54j l6Ct l6ll l63l l64l l65l l66l l67| . Given the current ambiguity on theoretical study of 
dark energy, in the observational probe of DE one often uses the parametrizations of EOS. 

Previously in the observational probes on oscillating features of dark energy EOS Ref. 68] has made some preliminary 
fittings to the pre-WMAP3 data and some relevant studies have been carried out later by Re fs.l53.|6 9l . In the present 
paper with the method dealing with the perturbations of Quintom developed in Refs. |l&l4!Al70ll71| . we aim to probe 
the time dependence of the dark energy EOS in light of WMAP3 and the combination with other tentative cosmological 
observations from SDSS and SNIa from the Riess "gold" sample or the SNLS observations. The background evolution 
and perturbations of Quintom can be identified with one normal quintessence and one phantom except for the phantom 
crossing point, where the natural matching condition is motivated by the case with a high-dimensional operator on the 
kinetic term or the two-field case and the individual sound speed for each field is assumed to be unity [Tsl 14^. Fol Fl | . 
In the present work we mainly focus on cases where the EOS is oscillating or with local bumps. By performing a 
global analysis with the Markov Chain Monte Carlo (MCMC) method, we find the current observations, in particular 
the WMAP3 + SDSS data combination, allow large oscillations of the EOS which can leave oscillating features on the 
(residual) Hubble diagram, and such oscillations are potentially detectable by future observations like SNAP. Local 
bumps of dark energy EOS can also leave imprints on CMB, LSS and SNIa. In cases when the bumps take place 
at low redshifts and the effective EOS is close to —1, CMB and LSS observations cannot give stringent constraints 
on such possibilities. However, geometrical observations like (future) SNIa can possibly detect such features. On 
the other hand when the local bumps take place at higher redshifts beyond the detectability of SNIa, future precise 
observations like Gamma-ray bursts and observations of 21 cm tomography, CMB and LSS may possibly detect such 
features. In particular, we find that bump-like dark energy EOS on high redshifts might be responsible for the features 
of WMAP on ranges / ~ 20 — 40, which is interesting and deserves addressing further. 

The remaining part of our paper is structured as follows: in Section II we describe the method and the data; in 
Section III we present our results on the determination of cosmological parameters with (WMAP3) Q, 0, S 0, , SNIa 
Sloan Digital Sky Survey 3-D power spectrum (SDSS-P(k)) by global fittings using the MCMC technique; 
discussions and conclusions are presented in the last section. 



In the parametrization of oscillating EOS one typically needs four parameters for the amplitude, center values, 
phase and the frequency. And in our analysis we have used 



The case with wq — — 1 and w\ = corresponds to the cosmological constant. 

The method we adopt is based on the publicly available Markov Chain Monte Carlo package CosmoMc l47ll7^ . which 
has been modified to allow for the inclusion of dark energy perturbations with EOS getting across — l|70j. Our most 




II. METHOD AND DATA 



w = wq + wi sin(w2 In a + W3) 



(1) 



1 For other solutions to this problem, see e.g. Refs. |41| - |45| and for a review on this issue see e.g. Ref. l46l . 
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general parameter space is 



p = (u! b ,uj c ,Qs,T,w Q , wi,w 2 , w 3 ,n s , log[10 




(2) 



where tut, = n^h 2 and to c — fl c h 2 are the physical baryon and cold dark matter densities relative to critical density, Os 
is the ratio (multiplied by 100) of the sound horizon to the angular diameter distance at decoupling, r is the optical 
depth, A s is defined as the amplitude of initial power spectrum and n s measures the spectral index. Assuming a flat 
Universe motivated by inflation and basing on the Bayesian analysis, we vary the above 10 parameters and fit to the 
observational data with the MCMC method. We take the weak priors as: r < 0.8, 0.5 < n s < 1 .5, —4 < wq < 1, — f < 
w\ < 1 0, < W2 < 20, — 7r/2 < w s < 7r/2, a cosmic age tophat prior as 10 Gyr< i <20 Gyr. The choice of priors on 
wq , u>i , u>2 , W3 have been set to allow for spread in all of the parameters simultaneously. Furthermore, we make use 
of the HST measurement of the Hubble parameter H — lOO/i km s^Mpc -1 [7^| by multiplying the likelihood by 
a Gaussian likelihood function centered around h — 0.72 and with a standard deviation a = 0.08. We impose a weak 
Gaussian prior on the baryon and density Slb/i 2 = 0.022 ± 0.002 (1 a) from Big Bang nucleosynthesis |7J. The bias 
factor of LSS has been used as a continuous parameter to give the minimum x 2 ■ 

In our calculations we have taken the total likelihood to be the products of the separate likelihoods of CMB, SNIa 
and LSS. Alternatively defining x 2 = — 21og£, we get 



In the computation of CMB we have included the three-year WMAP (WMAP3) data with the routine for computing 
the likelihood supplied by the WMAP team To be conservative but more robust, in the fittings to the 3D power 
spectrum of galaxies from the SDSS0| we have used the first 14 bins only, which are supposed to be well within the 
linear regime]]}. In the calculation of the likelihood from SNIa we have marginalized over the nuisance parameter |75|. 
The supernova data we use are the "gold" set of 157 SNIa published by Riess et al in 01 an d the 71 high redshift 
type la supernova discovered during the first year of the 5-year Supernova Legacy Survey (SNLS)|lg] resp ectively. In 
the fittings to SNLS we have used the additional 44 nearby SNIa, as also adopted by the SNLS group [l^- Also to be 
conservative but more robust, we did not try to combine SNLS with the Riess sample simultaneously for cosmological 
parameter constraints, namely in one case for SNIa fitting we use SNLS data only and in another case the Riess 
sample only. For each regular calculation, we run 6 independent chains comprising of 150,000-300,000 chain elements 
and spend thousands of CPU hours to calculate on a cluster. The average acceptance rate is about 40%. And for the 
convergence test typically we get the chains satisfy the Gelman and Rubin[7(j criteria where R-l<0.1. 

In our study for future perspectives on features of dark energy we have used the cosmic- variance |77) limited CMB 
TT spectrum up to I — 2000. For SNIa we have used SNAP[78l] simulations 2 and for LSS, we have adopted the 
LAMOST 80] simulations. In the remaining part of this paper the fiducial power law ACDM model adopted is as 
follows: 



where z r is the reionization redshift and the slightly different notations from previous Eq. Q are due to the difference 
in the CAMB f8lll82l | and CosmoMC^3lz3 default parameters. Such a fiducial model will be used to generate future 
CMB, SNIa and LSS data. In addition the illustrative figures will also be generated with such background parameters. 
Moreover in generating the illustrative figures on linear power spectrum of LSS we have fixed the bias factor to be 
unity. 

The projected satellite SNAP (Supernova / Acceleration Probe) would be a space based telescope with a one square 
degree field of view with 1 billion pixels. It aims to increase the discovery rate for SNIa to about 2,000 per year|78||. 
The simulated SNIa data distribution is taken from Refs. [831 l&l . l85 | . As for the error, we follow the ref. |83| which 
takes the magnitude dispersion 0.15 and the systematic error a sys (z) = 0.02 x z/1.7, and the whole error for each 
data is 



Xtotal — XCMB + XsNIa + XLSS 



(3) 



(u b , u c , h, z r , n s ,A s ) = (0.022, 0.12, 0.7, 12, 1, 2.3 x 10~ 9 ) , 



(4) 




(5) 



where n.; is the number of supernova in the i'th redshift bin. 



2 SNAP is one of the several candidate mission concepts for the Joint Dark Energy Mission (JDEM). Nowadays there have been many 
proposed dark energy surveys iTSll ■ 
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The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) project as one of the National Major 
Scientific Projects undertaken by the Chinese Academy of Science, aims to measure ~ 10 7 galaxies with mean redshift 
z ~ 0.2 80]. In the measurements of large scale matter power spectrum of galaxies there are generally two statistical 
errors: sample variance and shot noise. The uncertainty due to statistical effects, averaged over a radical bin Ak in 
Fourier space, is [8(| 

{ ^Lf = 2x ^l x 1 x(l + ^) 2 . (6) 
V P ' V 4vrfc 2 Afc V nP 1 1 ; 

The initial factor of 2 is due to the real property of the density field, V is the survey volume and n is the mean galaxy 
density. In our simulations for simplicity and to be conservative, we use only the linear matter power spectrum up to 
k ~ 0.15 h Mpc" 1 . 

For the cases with future cosmic variance limited CMB and LAMOST, we only show the error bars for illustrations, 
and a further analysis with fittings is not the aim of the present paper. For the studies on bump-like dark energy EOS, 
we also plot the illustrative results rather than make global fittings. And the parametrized EOS takes the following 
form: 

w = Wo + A(lna- A) 3 exp(-(lna- A) 4 /d) . (7) 

We should point out here that the parametrization in Eq. Q is a specific example only and there are certainly different 
parametrizations to illustrate the bump-like features in dark energy EOS and the resulting cosmological imprints might 
be different. 

III. RESULTS 

We start with the oscillating case. First of all in FigQ] we delineate the illustrative imprints of oscillating dark 
energy equation of state (EOS) on CMB (top left), LSS (top right) and on the Hubble diagram (lower panel). In the 
Hubble diagram the distance modulus \i is defined as the apparent magnitude m minus the absolute magnitude M: 



M = 51s <iWc) +25 ' (8) 



with dh being the luminosity distance: 

d L r dz' 



l + z J H(z') 



(9) 



In comparison the imprints by the ACDM cosmology have also been displayed. The main contribution of dark energy 
on CMB is on the geometrical angular diameter distance to the last scattering surface. This in turn determines the 
locations of CMB peaks. In some cases the effects on the large scale matter power spectrum of LSS are also significant. 
In Fig^we can find for cases where wo — —0.5, the effects on CMB turn out to be most eminently modulated. This 
is mainly due to the fact that in such cases as wq deviates significantly from —1 and w has been relatively matter- like, 
the contributions to large scale CMB are significant due to the Integrated Sachs- Wolfe (ISW) effects. An oscillating 
EOS of dark energy can leave somewhat similar imprints as oscillating primordial spectrum, as explored in Ref. |25| . 
although the difference is also noteworthy. We find interestingly that in cases wq = —1.5, the effects are not as large as 
wo = —0.5, which is in part due to the effects of dark energy perturbations and consistent with the previous analysis 
in Ref. [7(|. The effects of dark energy on CMB (and LSS) are mainly geometrical effects and can sometimes be 
understood from the formula on w e f / |87j| : 

_fdan(a)w(a) 

= fdan(a) • (10) 

On the other hand the Hubble diagram being displayed with the redshift, the effects of dark energy, in cases when 
DE dominates the Universe the EOS has oscillating features, can leave oscillations on the the Hubble diagram. 

For the next step we show the results from our global fittings on the oscillating EOS in Eq.lJIJ. In Table 1 we 
delineate the mean la constrains on the relevant cosmological parameters using different combination of WMAP, 
SNIa and SDSS. Shown together are those with maximum likelihood (dubbed ML) and those which give the most 
eminent oscillating effects in the 2a allowed regions (dubbed "Most" in the table). For simplicity the parameters 
related to bayon fractions and the primordial spectrum are not displayed, and the shown parameters are enough for 
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FIG. 1: Illustrative imprints of oscillating dark energy equation of state (EOS) on CMB (top left), LSS (top right) and on the 
Hubble diagram (lower panel) . 



the study on the Hubble diagrams below. Typically in the realizations of MCMC as the cosmological parameters are 
not exactly gaussian distributed, the center mean values are different from the best fit cases 3 . We find that although 
the best fit cases are given by oscillating EOS, a cosmological constant ( wq = —1, tui = ) is well within lcr for all 
the three data combinations. The accumulation of the observational data will help to break such a degeneracy. It is 
very interesting that wi is much better constrained by SNLS (together with the 44 low redshift SNIa) rather than by 
the Riess "gold" sample. On the other hand the background parameter Hq is better constrained by the Riess "gold" 
sample. One can also find that the frequency of oscillations, W2, is relatively the worst constrained parameter. We 
will turn to this in more details in the remaining part of the present paper. 

In Fig[21we delineate the one dimensional posterior constraints on the oscillating EOS in Eq.JI}, showing together 
some relevant background cosmological parameters. The black lines are constraints from combined analysis of WMAP3 
+ SDSS. The red lines are from WMAP3 + SDSS + Riess sample and the blue are from WMAP3 + SDSS + SNLS 
combined analysis. We can find that in some sense w% and are not well constrained by the current observations. As 
W3 represents the phase of oscillations and the range ( —tt/2, ir/2 ) is the largest 1-period limit (note we have allowed 
wi to be both positive and negative), our prior on u>2, though as large as given above, turn out to be somewhat too 
optimistic. 

The effects are more eminent in the resulting two dimensional contours. In Fig|2|we plot the corresponding two 
dimensional posterior constraints. The upper panel is on the constraints from the combined analysis of WMAP3 + 
SDSS. The lower left panel is from WMAP3 + SDSS + Riess sample and the lower right from WMAP3 + SDSS + 
SNLS combined analysis. 



3 More detailed discussions are available at http://cosmocoffee.info, 
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TABLE 1. Mean la constrains on cosmological parameters using different combination of WMAP, SNIa and SDSS. Shown 
together are those with maximum likelihood (ML) and those which give the most eminent oscillating effects within the 2a 

allowed regions (Most). 
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FIG. 2: One dimensional posterior constraints on the parametrized EOS: w — wo ± wi sin(w2 In a + U13) and on the relevant 
background cosmological parameters. The black solid lines are constraints from combined analysis of WMAP3 + SDSS. The 
red dash-dot lines are from WMAP3 ± SDSS ± Riess sample and the blue dashed lines are from WMAP3 ± SDSS ± SNLS 
combined analysis . 



It is understandable that in cases one or two parameters are not well constrained by the observations, changing 
the priors (in our case on W2) will inevitably affect the final results. Such a feature has also appeared in previous 
investigations on the oscillating primordial spectrum, regardless of fitting methods one uses such as grids |88j or 
MCMC [8f3.l90l| . For the grid case the problem is how small one should set the minimum grid, and there will be finer 
features on scales smaller than the minimum steps 4 . For our conventional cases with MCMC \18L l49l l7ll l91| to get the 
converged results, we test the convergence of the chains by Gelman and Rubin|7^] criteria and typically get R-l to be 
of order 0.01, which is more conservative than the recommended value R-l<0.1. However in the present case with an 
oscillating EOS, in the combinations with SNLS the maximum value of R-l is 0.04, for the combination with Riess 
"gold" sample, we have value: R-l ~ 0.1 and for the the case with WMAP3 ± SDSS we have the maximum R-l to 
be 0.25 instead, although for the last case we have run much longer time for the chains. This also shows that for the 
case with oscillations at least for WMAP3 ± SDSS our chains are not well converged. On the other hand this shows 
that for the current observations a large oscillation on the dark energy EOS is allowed. In Ref.[8£j the authors quoted 



4 For the grid case in the fittings to give the first year WMAP preferred running of the spectral index, in Ref. l24l in the case with an 
oscillating inflaton potential, the fractor-like structure is also present. Similar potentials can also lead to oscillations on the primordial 
spectrum, as previously investigated in Ref. l25l . And it is easily understood that in cases where the oscillations of the primordial 
spectrum have a large period, the resulting effects on CMB and LSS relevant scales will be running-like. Note besides inflation similar 
behavior also applies to the case of dark energy. 
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FIG. 3: Two dimensional posterior constraints on the parametrized EOS: w = wq + wi sin(w2 In a + IW3) at 68% and 95% C.L. 
The upper panel is on the constraints from combined analysis of WMAP3 + SDSS. The lower left panel is from WMAP3 + 
SDSS + Riess sample and the lower right from WMAP3 + SDSS + SNLS combined analysis . 

a value where R-l < 0.1 (their R-l is slig htl y different from ours) and the structures of our resulting two-dimensional 
figures are similar to those by Refs.[89L l9(| . These have implied that in the presence of oscillations typically the 
current data are not good enough to well break the degeneracy among the parameters. In our case although the 
parameters like n Sl wq, Qde, and Hq have already been well constrained, R-l is in some cases relatively large due 
to the parameters W2 and W3 . This has in turn led to the separate likelihood spaces in the three data combinations at 
68% and 95% C.L. Interestingly in the WMAP3 + SDSS + Riess sample the 1 a contours clearly separate into two 
different peaks in the lower three panels, resembling the contours in Ref. |8fl| . Such a behavior is less eminent in the 
WMAP3 + SDSS + SNLS sample, which nevertheless does exist in the W2 — W3 contour. We should point out that 
given the priors and the R-l value specified above, our results are robust. In Fig0]we delineate the resulting posterior 
lcr constraints on the low-redshift behavior of the oscillating EOS, with the three different data combinations. The 
red lines are given by the mean center values as shown in Table 1 and the blue dashes lines are the la allowed regions. 
The green dashed lines are the illustrative lcr explored regions by future SNAP[78j. Thus future SNIa observations 
like SNAP can help significantly to break the degeneracy, and in some sense detect such oscillating features of dark 
energy EOS. 

In Fig0 we delineate the corresponding imprints on the residue Hubble diagram (upper panels) and the Hubble 
diagram (lower panels). In the left panels, the lines dubbed "Osc" are given by the best fit values WMAP3 + SDSS, 
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FIG. 4: Resulting posterior constraints on the low-redshift behavior of the parametrized EOS: w = Wo + »i sin(w;2 In a + ws). 
The red lines are given by the mean center values as shown in Table 1 and the outside blue dashed lines are the la allowed 
regions by WMAP3 + SDSS (left), WMAP3 + SDSS + SNLS (center) and WMAP3 + SDSS + Riess sample (right). The 
inside green dashed lines are the illustrative la explored regions by future SNAP . 



the underlined "S" is by WMAP3 + SDSS + SNLS and "R" by WMAP3 + SDSS + Riess sample. The right panels 
show the corresponding cases where within the 2a allowed regions the oscillating effects are relatively eminent. Same 
as the above conventions in the upper panels the blue lines dubbed "SNAP" illustrate the detectability of future 
SNAP 78]. It is noteworthy that the combination with SNLS gives relatively the most stringent constraints. We can 
find that for the case with WMAP + SDSS only, a larger oscillation is allowed and there are some oscillating features 
on the (residue) Hubble diagram. And the effects of the "Most" case in Table 1 are more eminent, which could be 
possibly detected even by the CURRENTLY ONGOING SNIa observations. In the cases combined with Riess "gold" 
sample or with SNLS, the oscillating effects are less eminent. Nevertheless, oscillations are still present and in large 
areas of the parameter spaces, SNAP will be able to detect such features 5 . 

Now we turn to the studies on the local bump-like features of dark energy EOS. It is physically intuitive that dark 
energy EOS might not be exactly periodic, instead it could be semi-periodic or even with features on some specific 
redshifts. The parametrization adopted in Eq. Q) ( w — wq + ^4(lna — A) 3 exp(— (lna — A) 4 /<i) ) can accomodate 
constant EOS in cases where A = . A determines the locations of bumps, d determines the width and A determines 
the amplitudes of the bumps. One can easily understand that at the point lna = A we have w = wq and for lna < A 
(high redshifts) there will be a trough and for In a > A (low redshifts) there will be a peak (this happens only in cases 



5 In some sense in Table 1 for the "Most" case by WMAP + SDSS combination only, the relevant parameter space will get relatively 
stringently constrained by combinations with SNLS or with Riess "gold" sample. However given the possible discrepancies between 
SNLS and the Riess sample, which has been somewhat illustrated in Ref. Qj, currently ongoing SNIa observations can in some sense 
detect the (scmi)-oscillating features in dark energy EOS as allowed by WMAP + SDSS combination. 
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FIG. 5: Resulting imprints of the parametrized EOS: w — wo + wi sin(u)2 In a + W3) on the residue Hubble diagram (upper 
panels) and the Hubble diagram (lower panels). In the left panels, the lines dubbed "Osc" are given by the best fit values 
WMAP3 + SDSS, the underlined "S" is by WMAP3 + SDSS + SNLS and "R" by WMAP3 + SDSS + Riess sample. The 
right panels show the corresponding cases where within the 2a allowed regions the oscillating effects are relatively eminent. In 
the upper panels the blue lines dubbed "SNAP" illustrate the detectability of future SNAP . 



where A < 0, which is the region of our interest). And on very high redshifts where lno < A the second term will 
get damped exponentially and w approaches the value of wq. For simplicity and in the illustrative study here we fix 
wo = -1- 

Firstly we consider the cases where the bump takes place at low redshifts which can be detectable by (future) 
SNIa observations. In Fig|S] we delineate the imprints of bump-like dark energy EOS on the (residual) Hubble 
diagram, dark energy density fractions (left) and the corresponding effects on CMB and LSS (right). We have chosen 
A = 2 x 10 7 , A = -0.03 and d = 10~ 7 for the first illustration. From Fig© we can find that for such a specific 
choice of parameters the effects on CMB are indistinguishable from the corresponding ACDM model beyond cosmic 
variance, and nor can LSS tell one (bump-like) from the other (ACDM). It is understandable that in such a case as 
the effects of the peak and trough somewhat cancel out in the contributions to w e fj defined in Eq. 110t . the effects 
on CMB and LSS would then be almost indistinguishable from the ACDM case, and the results are consistent with 
those in Refs. [Tol. Wfy . In such a case, observations with geometric constraints on the redshift "tomography" will be of 
great importance to break such a degeneracy. The resulting (residual) Hubble diagrams are different from the ACDM 
model, and observations like SNAP can possibly detect such features. 

Although SNIa in some sense makes the only direct detection of dark energy, it is believed that due to some 
unavoidable effects like by dust the redshifts probed by SNIa cannot be too high, and for example in SNAP simulations 
one typically takes z < 1.7. In cases when the bump-like features take place at higher redshifts, we cannot expect 
to probe such features with even future SNAP observations. On the other hand as the current observations are 
consistent with a A-like dark energy, typically we cannot hope dark energy will be non-negligible in epochs far before 
the matter-radiation equality epoch, although such a possibility remains (for relevant studies see e.g. [92I l93| ) . On the 
other hand with our parametrized EOS given in Eq. J7J) we cannot expect dark energy component to be significant 
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FIG. 6: Imprints of bump-like dark energy EOS (denoted by "bulk model") on the (residual) Hubble diagram, dark energy 
density fractions (left) and the corresponding effects on CMB and LSS (right), shown together with the imprints by the ACDM 
model. The bump takes place at low redshifts which can be detectable by (future) SNIa observations. The effects on CMB are 
indistinguishable from the corresponding ACDM model beyond cosmic variance, and nor can LSS tell one (bump-like) from 
the other (ACDM). 



on very high redshifts, especially for cases where wq = — l 6 . Under such a circumstance we take the locations of 
bumps to be slightly larger than z = 2: in one case A = —1.5 and in another case A = —1.8. These correspond 
to z ~ 3.5 and 5.0 respectively where In a = A. For both of the cases we have fixed A = 6 x 10 5 and d = 10 -5 . 
In Fig[7|we show the imprints of bump-like dark energy EOS on CMB, LSS (right panels), on dark energy density 
fractions and on the residual Hubble diagram (left panels). The bumps take place at high redshifts which cannot be 
detectable by (future) SNIa observations. On the other hand, observations of Gamma-ray bursts (GRB) can probe 
relatively higher redshifts than SNIajgil HE HH > with the accumulations of GRB events as by SWIFT and better 
understandings on the systematics one is hopefully able to detect such features th rough geometric observations on 
the Hubble diagrams. Moreover observations like the 21 cm tomography |97ll98T. 193 . 1 100| are p otentially able to give 
almost a large number of independent samples compared with CMB and LSS observations 98J , features of bump- like 
dark energy EOS on intermediate redshifts are promisingly detected by observations of the 21 cm tomography. In 
the right panels of Fig[7| we find dramatically that bump-like features are possibly detectable by cosmic variance 
limited CMB (e.g. WMAP3HHH|l|i|) and LAMOST^. Part of the reason lies on the fact that with a different 
width and a higher redshift compared with the previous example in Fig|S| here the energy density fraction f2rj£;(a) 
increases and decreases rapidly and the magnitude of the peak value of Qde(o-) is somewhat lower, the asymmetric 
shapes of the peak and trough cannot compensate thoroughly to give w e ff close to —1. In fact one can find for our 
parametrization given in Eq. Q one will typically get w e ff > — 1 and in cases with larger —A the deviation of w e ff 
from -1 can be larger, which will also lead to some shifts on CMB peaks, as shown in Fig0 It is noteworthy to point 
out that both the first year WMAP and WMAP3 show some local glitches out of cosmic variance at I ~ 20 — 40. In 
our case of Fig0 there are also some relevant features on such scales and t his deserves fu rther investigations. While 
such features can be affected by different foreground analysis and so forth |l0ll Il02l ll03T| . it is theoretically possible 
that the features of WMAP TT on scales I ~ 20 — 40 are due to some bump-like features of dark energy EOS, or 
some semi-oscillations. Note here we have fixed the background parameters rather than performing a global analysis, 
and in cases we take into account the parameter degeneracies the detectability of the observations would typically get 



This is however, still possible given some fine tunings on the remaining parameters in Eq. J7J, which will not be explored by the current 
paper. 
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weaker. On the other hand we can expect that in the global analysis when all of the parameters (including the bias 
factor) can vary, the observations can give constraints or detect the signatures where the amplitudes of bumps are 
even higher than the examples listed above. 




FIG. 7: Imprints of bump-like dark energy EOS on CMB, LSS (right panels), on dark energy density fractions and on the 
residual Hubble diagram (left panels). The bumps take place at high redshifts which cannot be detectable by (future) SNIa 
observations. On the other hand, cosmic variance limited CMB (e.g. WMAP3) and LAMOST|8(j| are possible to detect such 
features. Note here we have not performed a global analysis and did not take into account the parameter degeneracies . 



IV. DISCUSSION AND CONCLUSION 



We should point that in all our three data combinations, a large area of oscillating Quintom where the EOS of 
dark energy gets across —1 during time evolutions are allowed by the current observations , w hich is different from 
the previous case in Ref.[5(| and from the case in Ref. (53J- Moreover contrary to Refs.jiO, E3, oscillating Quintom 
displays the distinctive feature that due to the average effects on the regions of EOS where w > — 1 and w < — 1, the 
quantity iv e ff can be very close to —1. Oscillating Quintom- like dark energy, where the unification of two epochs of 
accelerated expansions can be realized and the coincidence problem can be solved jB^, proves to be a good fit to current 
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cosmology. We should stress again that resembling the oscillating primordial spectrum case[8!j, the not good enough 
convergence of oscillating dark energy EOS is inevitable in the light of the current observations. In the WMAP + SDSS 
allowed parameter space there are areas with significant oscillations, which implies that the CURRENTLY ONGOING 
SNIa observations may detect such oscillating or semi-oscillating features on the (residual) Hubble diagram. Moreover 
while oscillating features on CMB might be explained by both oscillating primordial spectrum and oscillations in dark 
energy EOS, the oscillating features on the Hubble diagram cannot be due to features in the primordial spectrum, 
but by oscillating EOS. Local bumps of dark energy EOS may leave some distinctive imprints on CMB, LSS and 
SNIa. The bumps are potentially detectable by geometrical observations like SNIa and GRB. Future observations of 
21 cm tomography open a very promising window to detect/exclude such features. LSS measurements like LAMOST 
and cosmic variance limited CMB may also detect such features. In particular, bump-like dark energy EOS on high 
redshifts might be responsible for the features of WMAP on ranges I ~ 20 — 40, which is interesting and deserves 
addressing further. 

Given the fact that currently we know relatively very little on theoretical aspects of dark energy, observing dynamical 
dark energy is currently the most important aspect of dark energy study. On the observational probes typically one 
needs some parametrizations on the form of dark energ y EOS. While using some specific forms of pa rametrizations 
one often takes the risks of getting biased ( see e.g. |104| ). investigations towards unbiased probe of DE |l05lfl06lll07| 
in light of all the available cosmological observations still need f urthe r developments. On the other hand some of 
the parametrizations are in some sense well motivated (see e.g. |104| |'). and in fittings starting with parametrized 
EOS one typically has larger v (number of data minus the number of parameters) and hence gets better constraints 
on dark energy. With the accumulations of observational data poor parametrizations will get ruled out and in this 
sense parametrized study of dark energy provides a complementary study of the non-parametric probes. As fro m an 
quintessence-like or phantom- like EOS which do not get across —1 one can reconstruct the dark energy potential|lQi 
any parametrizations of non-Quintom-like EOS correspond to some specific forms of quintessence/phantom potentials 7 
and hence such parametrizations are somewhat well motivated. Bump- like EOS which do not get across —1 also 
correspond to some specific DE potentials which can be straightforwardly worked out. Similarly for (semi-)oscillating 
EOS which do not get across —1 one can work out the corresponding DE potentials. In cases with bumps described 
in Eq. (JJJ and with oscillations in Eq. (JIJ , one cannot simply reconstruct the potentials of DE due to the distinctive 
nature of Quintom[il|7(3. However this r emains possible for exam ple the models with high derivatives [6^, |(37| and 
in the framework with modified gravity [H EM El3 EH EH EH • 

Acknowledgments: We acknowledge the use of the Legacy Archive for Microwave Background Data Analysis 
(LAMBDA). Support for LAMBDA is provided by the NASA Office of Space Science. We have performed our 
numerical analysis on t he Shangh ai Supercomputer Center (SSC). We used a modified version of C AMB |8lL l82| which 
is based on CMBFAST pll 115]. We are grateful to Hiranya Pciris, Yunsong Piao and Lifan Wang for discussions 
related to this project. We thank Xuelei Chen, Yaoquan Chu and Long-long Feng for discussions on LAMOST. We 
thank Sarah Bridle, Peihong Gu, Steen Hannestad, Antony Lewis, Mingzhe Li, Yongzhong Xu, Jun'ichi Yokoyama, 
Max Tegmark and Penjic Zhang for helpful discussions. B. F. would like to thank the hospitalities of IHEP during his 
visit to Beijing. This work is supported in part by National Natural Science Foundation of China under Grant Nos. 
90303004, 10533010 and 19925523 and by Ministry of Science and Technology of China under Grant No. NKBRSF 
G19990754. B. F. is supported by the JSPS fellowship program. 



D. N. Spergel et al.,|arXiv:ast ro-ph/0603449 
L. Page et al, arXiv:astro-ph/0603450 
G. Hinshaw et al, arXiv:astro-ph/0603451 

N. Jarosik et al, arXiv :astro-ph /06034 52 | 

Available from http://lambda.gsfc. nasa.gov/product/map/current/]- 

M. Tegmark et al. [SDSS Collaboration], Astrophys. J. 606, 702 (2004) |arXiv:astro^ph/ 0310725 . 
M. Tegmark et al. [SDSS Collaboration], Phys. Rev. D 69, 103501 (2004) |arXiv:astro-pfi/0310723| . 
S. Weinberg, Rev. Mod. Phys. 61, 1 (1989). 



I. Zlatev, L. M. Wang and P. J. Steinhardt, Phys. Rev. Lett. 82, 896 (1999) arXiv:astro-ph/9807002 . 
R. D. Peccei, J. Sola and C. Wetterich, Phys. Lett. B 195, 183 (1987); C. Wetterich, Nucl. Phys. B 302, 668 (1988); 

B. Ratra and P. J. E. Peebles, Phys. Rev. D 37, 3406 (1988). 

[11] A. G. Riess et al. [Supernova Search Team Collaboration], Astron. J. 116, 1009 (1998) [arXiv:ast ro-ph/9805201 . 



In cases where w > 1 they correspond to quintessence with negative potentials. 



13 



[12] S. Perlmutter et al . [Supernova Cosmology Project Collaboration], Astrophys. J. 517, 565 (1999) 
|arXiv:astro-ph/9812133] . 

J. L. Tonry et al. [Supernova Search Team Collaboration], Astrophys. J. 594, 1 (2003) arXiv:astro-ph/0305008 . 
A. G. Riess et al. (Supernova Search Team Collaboration), Astrophys. J. 607, 665 (2004) arXiv:astro-ph/0402512 . 

A. Clocchiatti et al. (the High Z SN Search Collaboration), astro-ph/0510155 
P. Astier et al, Astron. Astrophys. 447, 31 (2006) arXiv:astro-ph/05 10447 . 
U. Seljak, A. Slosar and P. McDonald, arXiv:astro-ph/0604335 

G. B. Zhao, J. Q. Xia, B. F eng and X. Zhang, |arXiv:astro-p h/0603621 
Y. Wang and P. Mukherjee, arXiv:astro-ph/0604051 

C. L. Bennett et al, Astrophys. J. Suppl. 148, 1 (2003) [arXiv:astro-ph/0302207| . 

D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 148, 175 (2003) arXiv:astro-ph/0302209 . 

H. V. Peiris et al, Astrophys. J. Suppl. 148, 213 (2003) |arXiv:astro -ph /0302225 . 

e.g. C. R. Contaldi, M. Peloso, L. Kofman and A. Lind^TJCAP 0307, 002 (2003 ); B. Feng and X . Zhang, Phys. 
Lett. B 570, 145 (2003); Q. G. Huang and M. Li, JCAP 0311, 001 (2003) |arXiv:astro- ph /0308458 ; M. Yamaguchi 
and J. Yokoyama, Phys. Rev. D 70 (2004) 023513; Y. S. Piao, B. Feng and X. m. Zhang, Phys. Rev. D 69, 103520 
(2004 ) |arXiv:h ep-th/0310206 ; Y. S. Piao, S. Tsujikawa and X. m. Zhang, Class. Quant. Grav. 21, 4455 (2004) 
|arXiv:hep-th/0312139|; Y. S. Piao, Phy s. Rev. D 71, 087301 (2005) |arXiv:astro^ph /0502343 ; B. A. Bassett, S. Tsujikawa 
and D. Wands, |arXiv:astr o-ph/0507632 

B. Feng, M. z. Li, R. J. Zhang and X. m. Zhang, Phys. Rev. D 68, 103511 (2003) arXiv:astro-ph/0302479 . 
X. Wang, B. Feng, M. Li, X. L. Chen and X. Zhang, Int. J. Mod Phys. D 14, 1347 (2005 ) |arXiv:astro-ph/0209242| . 
T. Moroi and T. Takahashi, Phys. Rev. Lett. 92, 091301 (2004) arXiv:astro-ph/0308208 . 
B. Feng, X. Gong and X. Wang, Mod. Phys. Lett. A 19, 2377 (2 004) |arXiv:astro-ph/0301111| . 
M. Viel, M. G. Haehnelt and A. Lewis, arXiv:astro-ph/0604310 
L. A. Kofman, A. D. Linde, A. A. Starobinsky, Phys. Lett. B 157, 361 (1985). 

A. A. Starobinsky, JETP Lett. 55, 489 (1992) [Pisma Zh. Eksp. Teor. Fiz. 55, 477 (1992)]. 

J. A. Adams, G. G. Ross and S. Sarkar, Nucl. Phys. B 503, 405 (1997) |arXiv:hep-ph/9704286[. 

J. Lesgourgues, D. Polarski and A. A. Starobinsky, Mon. Not. Roy. Astron. Soc. 297, 769 (19 98) |arXiv:astro-ph/9711139"] . 
D. J. H. Chung, E. W. Kolb, A. Riotto and I. I. Tkachev, Phys. Rev. D 62, 043508 (2000) arXiv:hep-ph/9910437 . 
L. M. Wang and M. Kamionkowski, Phys. Rev. D 61, 063504 (2000) arXiv:astro-ph/9907431 . 

J. Lesgourgues, Nucl. Phys. B 582, 593 (2000) arXiv:hep-ph/9911447 . 

L. M. Griffiths, J. Silk and S. Zaroubi, Mon. Not. Roy. Astron. Soc. 324, 712 (2 001) |arXiv:astro-ph/0010571| . 
S. Hannestad, S. H. Hansen and F. L. Villante, Astropart. Phys. 16, 137 (2001) arXiv:astro-ph/0012009 . 
J. Barriga, E. Gaztan aga, M. G. Santos and S. Sarkar, Mon. Not. Roy. Astron. Soc. 324, 977 (2001) 

arXiv:astro-ph/0011398 . 

M. Gramann and G. Hutsi, Mon. Not. Roy. Astron. Soc. 327, 538 (2001) arXiv:astro-ph/0102466 . 

M. Kamionkowski and A. R. Liddle, Phys. Rev. Lett. 84, 4525 (2000) arXiv:astro-ph/9911103 ; A. R. Zentner and 

J. S. Bullock, Phys. Rev. D 66, 043003 (2002) [arXiv:astro-ph/0205216l . 

D. N. Spergel and P. J. Steinhardt, Phys. Rev. Lett. 84, 3760 (2000) arXiv:astro-ph/9909386 . 

B. D. Wandelt, et al. Proceedings of Dark Matter 2000, |arXiv:astro-pli/0006344| 

M. Kaplinghat, L. Knox and M. S. Turner, Phys. Rev. Lett. 85, 3335 (2000) " |arXiv:astrc-p h/0005210]. 

W. B. Lin, D. H. Huang, X. Zhang and R. H. Brandenberger, Phys. Rev. Lett. 86, 954 (2001 ) |arXn^ :astro-ph/0009003 1. 
P. Bode, J. P. Ostriker and N. Turok, Astrophys. J. 556, 93 (2001) arXiv:astro-ph/0010389 . 
A. Tasitsiomi, Int. J. Mod. Phys. D 12, 1157 (2003) arXiv:astro-ph/0205464 . 

A. Lewis and S. Bridle, Phys. Rev. D 66, 103511 (2002) |arXiv:astro-ph/0205436|. 

P. J. E. Peebles and A. Vilenkin, Phys. Rev. D 59, 06350 5"TT999) |arXiv:astro-ph /9810509 . 
J. Q. Xia, G. B. Zhao, B. Feng and X. Zhang, arXiv:astro-ph/0603393 

S. Dodelson, M. Kaplinghat and E. Stewart, Phys. Rev. Lett. 85, 5276 ( 2000). 

B. Feng, M. Li, Y. S. Piao and X. Zhang, Phys. Lett. B 634, 10 1 (2006) |arXiv:astro-ph/0407432| . 
G. Barenboim and J. D. Lykken, Phys. Lett. B 633, 453 (2006) |arXiv:astrcTp h/0504090 . 

For an interesting relevant study see G. Barenboim, O. Mena and C. Quigg, Phys. Rev. D 71, 063533 (2005) 
arXiv:astro-ph/0412010 . 

B. Feng, X. L. Wang and X. M. Zhang, Phys. Lett. B 607, 35 (2005) arXiv:astro-ph/0404224 . 
S. Weinberg, Phys. Rev. Lett. 59, 2607 (1987). 

G. R. Dvali, G. Gabadadze and M. Porrati, P hys. Let t. B 485, 208 (2000) [arXiv:h ep-th/0005016 . 
K. Koyama, Phys. Rev. D 72, 123511 (2005) arXiv:hep-th/0503191 . 

For an interesting relevant study see J. Yokoyama, Phys. Rev . Lett. 88, 151302 (2002) [arXiv:hep-th/0110137l . 
For a recent reivew see A. Vilenkin, |arXiv:astro-ph/0605242| 

For a review see E. J. Copeland, M. Sami and S. Tsujikawa, arXiv:hep-th/0603057 

e.g. H. Wei and R.-G. Cai, arXiv:hep-th/0501160 R.-G. CaiHl.-S. Zh ang and A . Wang, |arXiv:hep-th/0 505186 A. A. 
Andrianov, F. Cannata and A. Y. Kamenshchik, arXiv:gr-qc/0505087 X. Zhang, arXiv:astro-ph/0504586 Q. Guo and 



Vernov, arXiv:astro-ph/0507067 


C. G. Huang and H. Y. Guo, arXiv:astro-ph/0508171 W. Zhao, 


arXiv:astro-ph/0604460 


J. Grande, J. Sola and H. Stefancic, 


arXiv:gr-qc/0604057 


L. P. Chimento and R 


Lazkoz, 


arXiv:astro-ph/0604090 


F. Cannata and A. Y. Kamenshchik, arXiv:gr-qc/0603129 


R. Lazkoz and G. Leon, arXiv:astro-ph/0602590 


B. Feng, 



14 



arXiv:astro-ph/0602156 H. Stefancic, arXiv:astro-ph/0511316 L. Perivolaropoulos, arXiv:astro-ph/0601014 
[62] T. Chiba, T. Okabe and M. Yamaguchi, Phys. Rev. D 62 (2000) 023511 |arXiv:astro-ph/9912463| ; C. Armendariz-Picon, 
V. F. Mukhanov and P. J. Steinhardt, Phys. Rev. L ett. 85, 4438 (2000) |arXiv:astro -ph/0004134 . 

R. R. Caldwell, Phys. Lett. B 545, 23 (2002) arXiv:astro-ph/9908168 . 

Z. K. Guo, Y. S. Piao, X. M. Zhang and Y. Z. Zhang, Phys. Lett. B 608, 177 (2005) arXiv:astro-ph/0410654 . 

M. z. Li, B. Feng and X. m. Zhang, JCAP 0512, 002 (2005) [arXiv:hep-ph/0503268| . 

X. F. Zhang, H. Li, Y. S Piao and X. M. Zhang, Mod. Phys. Lett. A 21, 231 (2006) |arXiv:astro-'pl7/ 0501652 . 

X. F. Zhang and T. Qiu, arXiv:astro-ph/0603824 

J. Q. Xia, B. Feng and X. M. Zhang, Mod. Phys. Lett. A 20, 2409 (2005) arXiv:astro-ph/0411501 . 
E. V. Linder, Astropart. Phys. 25, 167 (2006) [arXiv:astro-ph/0511415l . 

G. B. Zhao, J. Q. Xia, M. Li, B. Feng and X. Zhang, Phys. Rev. D 72, 123515 (2005) arXiv:astro-ph/0507482 . 
J. Q. Xia, G. B. Zhao, B. Feng, H. Li and X. Zhang, Phys. Rev. D 73, 063521 (2006) |arXiv:astro-ph/051 1625J . 
Available from http://cosmologist.info. 

W. L. Freedman et al, Astrophys. J. 553, 47 (2001) |arXiv:astro-ph/0012376"1. 

S. Buries, K. M. Nollett and M. S. Turner, Astrophys. J. 552. I.I (21101) |arXiv:astro-ph/0010171| . 
For deta ils see e.g. E . Di Pietro and J. F. Claeskens, Mon. Not. Roy. Astron. Soc. 341, 1299 (2003), 
arXiv:astro-ph/0207332 . 

A. Gelman and D. Rubin, Statistical Science 7, 457 (1992). 
J. R. Bond and G. Efstathiou, Mon. Not. Roy. Astron. Soc. 226, 655 (1987). 
Available at http://snap.lbl.g ov] . 

See e.g. A' Crotts ei aC, |arXiv:astro-ph/0507043| T. Abbott et al. [Dark Energy Survey Collaboration], 
ar Xiv:astro-ph/05 10346 ] 

Available at http://www.lamost.org/ . 

A. Lewis, A. Challinor a nd A Lasenby, Astrophys. J. 538, 473 (2000) arXiv:astro-ph/9911177 . 

Available at http://camb.info . 

A. G. Kim, E. V. Linder, R. Miquel and N. Mostek, Mon. Not. Roy. Astron. Soc. 347, 909 (2004) arXiv:astro-ph/0304509 . 

C. Yeche, A. Ealet, A. Refregier, C. Tao, A. Tilquin, J. M. Virey and D. Yvon, arXiv:astro-ph/0507170 

H. Li, B. Feng, J. Q. Xia and X. Zhang, Phys. Rev. D 73, 103503 (2006) arXiv:astro-ph/0509272 . 

H. A. Feldman, N. Kaiser and J. A. Peacock, Astrophys. J. 426, 23 (1994) |arXiv:astro-ph/9304022|. 

e.g. L. M. Wang, R. R. Caldwell, J. P. Ostriker and P. J. Steinhardt, Astrophys. J. 530, 17 (2000) arXiv:astro-ph/9901388 . 
O. Elgaroy and S. Hannestad, Phys. Rev. D 68, 123513 (2003) |arXiv:astro- ph/03070li"1. 
T. Okamoto and E. A. Lim, Phys. Rev. D 69, 083519 (2004) |arXiv:astro-ph/0312284T~^ 
R. Easther, W. H. Kinney and H. Peiris, JCAP 0505, 00 9 (2005) |arXiv:astro-ph/0412613l . 

B. Feng, M. Li, J. Q. Xia, X. Chen and X. Zhang, arXiv:astro-ph/0601095 
M. Malquarti and A. R. Liddle, Phys. Rev. D 66, 023524 (2002) arXiv:astro-ph/0203232 . 
W. Wang and B. Feng, Chin. J. Astron. Astrophy s. 3, 105 (2003) |arXiv :astro-ph/0508139 . 

K. Takahashi, M. Oguri, K. Kotake and H. Ohno, arXiv:astro-ph/0305260 

Z. G. Dai, E. W. Liang and D. Xu, Astrophys. J. 612 , L101 (2004) |arXiv:astro-ph /0407497 . 

D. Hooper and S. Dodelson, arXiv:astro-ph/0512232 

G. B. Field, Astrophys. J. 129, 536 (1959). 

A. Loeb and M. Zaldarriaga, Phys. Rev. Lett. 92, 211301 (2004) |arXiv:astro-ph/0312134| . 
X. L. Chen and J. Miralda-Escude, Astrophys. J. 602, 1 (200 4) |arXiv:astro -ph/0303395r 
e.g. U. L. Pen, X. P. Wu and J. Peterson, arXiv:astro-ph/0404083 

See also A. Slosar, U. Seljak and A. Makarov, "Exact likelihood evaluations and foreground marginalization in low 
resolution WMAP data," Phys. Rev. D 69, 123003 (2004) arXiv:astro-ph/0403073 . 

A. de Oliveira-Costa and M. Tegmark, Phys. Rev. D 74, 023005 (2006) larXiv:ast ro-ph/0603369l. 

H. K. Eriksen et al., "A re-analysis of the three-year WMAP temperature power spectrum and arXiv:astro-ph/0606088 

B. A. Bassett, P. S. Corasaniti and M. Kunz, Astrophys. J. 617, LI (2004) arXiv:astro-ph/0407364 . 
D. Huterer and G. Starkman, Phys. Rev. Lett. 90, 031301 (2003) arXiv:astro-ph/0207517 . 
Y. Wang and M. Tegmark, Phys. Rev. D 71, 103513 (2005) |arXiv:astro-ph/0501351|. 
J. Simon, L. Verde and R. Jimenez, Phys. Rev. D 71, 123001 (2005) |arXiv:astro-p h/0412269 . 
A. Lewis, CAMB notes. 

M. Morikawa, Astrophys. J. 362, L 37 (1990); Astr ophys. J. 3 69, 20 (1991). 
V. Sahni and Y. Shtanov, JCAP 0311, 014 (2003) arXiv:astro-ph/0202346 . 
L. Perivolaropoulos, Phys. Rev. D 67, 123516 (2003) |arXiv:hep-ph/0301237l . 
H. Stefancic, Eur. Phys. J. C 36, 523 (2004) arXiv:astro-ph/0312484 . 

L. Perivolaropoulos, JCAP 0510, 001 (2005) |arXiv:astro-ph/050458'2|. 

U. Seljak and M. Zaldarriaga, Astrophys. J. 469, 437 (1996) arXiv:astro-ph/9603033 . 
Available at http://cmbfast.org/ . 



